Controllable spin-current blockade in a Hubbard chain 
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We investigate the spin/charge transport in a one-dimensional strongly correlated system by using 
the adaptive time-dependent density-matrix renormalization group method. The model we consider 
is a non-half-filled Hubbard chain with a bond of controllable spin-dependent electron hoppings, 
which is found to cause a blockade of spin current with little influence on charge current. We have 
considered (1) the spread of a wave packet of both spin and charge in the Hubbard chain and (2) 
the spin and charge currents induced by a spin-dependent voltage bias that is applied to the ideal 
leads attached at the ends of this Hubbard chain, ft is found that the spin-charge separation plays 
a crucial role in the spin-current blockade, and one may utilize this phenomenon to observe the 
spin-charge separation directly. 

PACS numbers: 72.25.-b, 71.10.Pm, 71.10.Pd, 73.23.Hk 



In one-dimensional (Id) strongly correlated systems, 
an essential phenomenon is the spin-charge separation 
(SCS), [H which is believed to play a central role in Id 
transport Q. To study transport problems in Id sys- 
tems, the characteristics of SCS must be taken into ac- 
count. However, due to the limitations of existing meth- 
ods, the discussion in the past has been limited to a few 
simple cases: an example is the significant work by Kane 
and Fisher [3[ on scaling properties of tunneling through 
a spin-symmetric point impurity in a fermion system. 
Hence more powerful methods are needed to compute 
transport properties beyond scaling and to treat a general 
interacting Hamiltonian. Furthermore, spin-dependent 
transport problems have attracted increasing interest in 
the last two decades: many proposed spintronics de- 
vices 0, H[ require the manipulation of spin currents, 
which are naturally decoupled from charge currents in Id 
systems. On the other hand, a number of experimental 
works have sought to observe the phenomenon. Rapid 
progress in ultracold atomic gas experiments Q makes it 
possible to see SCS in a new context. [8| 

Recently, the adaptive time-dependent density-matrix 
renormalization group 

(t-DMRG)HJli 

was developed 

by combining the DMRG met hod 1 1 with quantum 
information concepts. The key idea of this method 
is to break u p t he evolution operator with Trotter 
decomposition, 12] then apply it to the states within a 
DMRG configuration. At the same time, other real- 
time evolution methods within DMRG[H, Q [H, [H 
were also proposed. By use of these numerical methods, 
there have been investigations on transport properties 
in Id strongly-correlated or impurity systems, including 
spin- 1/2 chains. [l7l B ose- Hubbard model. (lq and con- 
ductance analysis. 19] The dynamical problems with im- 
purities were also studied widely using static DMRG 
method embedding with persistent current [i^, [2l| and 
functional renormalization group. [2lL 22 ] An interesting 



result that partly motivates our study was the study of 
SCS by Kollath, Schollwoeck, and Zwerger.jiij] In con- 
ventional treatment with the bosonization method, only 
low energy excitations were considered. The above study 
goes beyond the low energy excitation spectrum by con- 
sidering the evolution of a "big" (multiparticle) wave 
packet that shows the SCS phenomenon. [231] 

In this Letter, we propose to consider a non-half- filled 
Hubbard chain in which one special bond has control- 
lable spin-dependent electronic hoppings, motivated by 
the development of optical lattices of ultracold atoms 
in which all hoppings can be controlled. By using the 
adaptive t-DMRG method, we simulate the spread of 
a wave packet as well as the spin and charge currents 
under a spin-dependent voltage bias. We find that the 
spin-current blockade can be realized by adjusting the 
spin-dependent hopping on that special bond while the 
charge current has not been affected. A possible appli- 
cation of this phenomenon is discussed. 

The system we are considering is described by the fol- 
lowing Hamiltonian, 



+ h.c.) + [/2^n a n u , (1) 



where cj CT (ci )0 -) creates (annihilates) an electron with 

spin <j(=T ) I) on the i-th site, n,i y(T (= c\ a Ci. a ) is the 
corresponding electron number operator, the electron 
hopping constants tf i+1 = to on all bonds but a spe- 
cial bond (i = l s ), where tf i+1 (= t a ) is an adjustable 
spin-dependent quantity, which introduces a local mag- 
netic moment, like a spin-dependent Anderson impurity. 
Without the special bond, the Hamiltonian is nothing but 
the usual Hubbard model with t being the electron hop- 
ping constant between nearest neighbor sites and U(> 0) 
the on-site repulsive Coulomb interaction. Without loss 
of generality, we will always keep t-\ = to while t± is ad- 
justed from to to 0. We present a proposal in Fig. 1 on 
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FIG. 1: The realization of a Hubbard chain with a special 
bond of controllable spin-dependent electronic hopping in an 
optical lattice. The rotation of the left spin-down atoms 
around the spin-up atom axis reduces the spin-down atom 
hopping at the special bond, i.e., the t± in Eq. (1) decreases 
as the angle 6 increases. 



the realization of a Hubbard chain with a special bond 
of controllable spin-dependent electronic hopping in an 
optical lattice. In the optical lattice, the spin-up and 
-down atoms could be trapped in a pair of parallel pe- 
riod potential wells created by interfering laser beams. [Jj 
Hence, one might change the laser beam to rotate half 
of spin-down atoms around the spin-up atom axis by an 
angle 9, which reduces the spin-down atom hopping at 
the special bond. It is clear that of this special bond 
in Eq. (1) could be easily adjusted by changing 9. 

First, we consider the spread of a wave packet as done 
by Kollath et al. [23j, but in a Hubbard chain with a 
bond of controllable spin-dependent electron hoppings, 
described by Eq.(l). The wave packet is induced by the 
following spin-dependent Gaussian potential 



Hp = -P^exp 



(» - y 2 



21/ 



H,Ti 



(2) 



which acts only on the spin up electrons, so that it car- 
ries both of spin and charge. Then the potential will be 
switched off and the wave packet will be spread according 
to the time-dependent many-body Schrodinger equation. 
Clearly P determines the potential strength, l p the cen- 
ter of the induced wave packet and Id its width. A small 
value of Id (e.g., <1) corresponds to a local potential at a 
single site as suggested by Anderson, the case we will 
consider here. 

In the following, we apply the t-DMRG method with 
second order Trotter decomposition to simulate the dy- 
namical evolution of this wave packet. The number of 
sites (L) is taken to be 176 while the number of electrons 
(N e ) is 116; the corresponding filling factor (N e /(2L)) is 
about 1/3. The time step is taken as 0.04 (in unit of h/to) 
and the number of kept DMRG states (M) is chosen to be 
large enough to ensure the error being less than O(10 -3 ). 
The spin and charge currents, the quantities we focus on 
in this calculation, are defined as J s {j) = J](j) — ■Tj.O) 
and J c (j) = J-\(j) + J[{j), separately, where 



(3) 



In Fig. [21 we show the spin and charge current at vari- 
ous times for = . The initial potential locally on the 
spin up generates simultaneously spin and charge wave 
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FIG. 2: The spin (solid line) and charge (dashed line) current 
in the spread of a wave packet in a Hubbard chain with a 
special bond with U = Sto, l s = 83 and P — 0.8to,ld = 1, lp = 
62. The location of this special bond is indicated by a double 
dashed (blue) line. Positive values imply the current flows 
to the right while negative values the current flows to the 
left. The dashed (red) arrow indicates the charge transmitted 
through the special bond while the solid (black) one the spin 
reflected from the special bond. 



packet, and then the wave packets will split into two parts 
and propagate to opposite directions respectively. [23j 
Since the potential strength P we took is not large, the 
spin and charge density of this wave packet is very small, 
but the currents defined in Eq.(3) show the propagation 
of this wave packet very clear. In Fig. HJa), we show the 
current distribution at a time before the spin and charge 
reach the special bond. The split of the two peaks indi- 
cates the different speeds of spin and charge excitations, 
so that the spin-charge separation is observed clearly. 
Next we'll see clearly from the figure, the charge cur- 
rent goes through the special bond almost freely, while 
the spin current is blocked by the special bond. A spin 
current reflecting at the special bond is shown in the fig- 
ure by its value being changed from positive to negative. 
The change of charge currents from negative to positive 
indicates the charge reflection at the left end since we use 
an open boundary condition for the chain. 

We argue that, in the following, the spin-current block- 
ade we observe here only happens in a strongly corre- 
lated svstem. It is different from the "spin blockade" 
effect, [24| which is merely spin-related Coulomb block- 
ade. To understand the phenomenon we observed above, 
we consider the large U limit of the Hubbard-like model 
in Eq. (1), which leads to the so-called t — J model 
H t -j = Ht + Hj, where H t is the hopping term in Eq. (1) 
and 

H.J = -77 ^ ^jAsCjsn^njicl^Cks, (4) 



ijkss' 



with t? 



25] This model works 



on the space that has projected out all configurations 
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FIG. 3: A spin-exchange process that is necessary for the spin 
current in the strong interaction limit. 



with at least one doubly occupied site for a less-than- 
half filling system and is responsible for the low-energy 
excitations of the model in Eq. (1). It is well known 
that the hopping term H t is responsible for the charge 
excitations while Hj controls the spin excitations, and 
that for the usual Hubbard model, Hj corresponds to a 
Heisenberg spin chain. Writing 



Si ' Si — Si Si 



2 i^i Sj 



S~ S ' - 



(5) 



the second term is responsible for the spin-exchange pro- 
cess that is necessary for an S z spin current. Based on 
Hj in Eq.(4) from the large-C/ expansion, we show the 
spin-exchange process in Fig. 3, in which there is a virtual 
intermediate state, so electronic hoppings for both spins 
are necessary for the process via a virtual state. Then it 
is clear that spin current is blocked at the special bond 
when t[ = since only spin-up electrons can hop across 
the special bond. 

Now, we come to the calculation of spin and charge 
currents under a spin-dependent voltage bias, for which 
we consider the system of a Hubbard chain with a special 
bond is attached with two ideal leads at its two end, that 
is described by the following Hamiltonian, 



H = H, 



H 



lead 



(6) 



where H$ has been given in Eq.(l) and the summation 
over site index runs from 1 to L with L being the size of 
the Hubbard chain, Hi ea d = Hi, + Hr and 



H n — —tr 



(7) 



with a — L or R and the sign ± taking — (+) for a 
L (R), and 



Hint — ~ to 



^( c (L C l,ff + A^°L+l,a + h-C). (8) 



The spin-dependent voltage bias V a = — Vj 1 applied 
at the two leads turns on at t — 0, so that currents appear 
gradually at the same time and reach constant values 
finally. In the calculation, we take Vi = and V\ = O.lin, 
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FIG. 4: The spin and charge currents at various times for dif- 
ferent numbers of kept DMRG states (M) or Hubbard chain 
lengths (L). (a) Both incoming and outgoing spin/charge cur- 
rents are shown for {7 = and L = 40. (b) Only the outgoing 
currents are shown for U = 8to and L = 40. (c) The outgo- 
ing currents for different chain lengths L, and U = 8io and 
M = 160. 



without loss of generality, to induce both spin and charge 
currents through the Hubbard chain. 

Now we show in Figure 4 the spin and charge currents 
at various times for different Hubbard chain lengths (L) 
or the numbers of kept DMRG states (M) in the calcula- 
tions induced by the spin-dependent voltage bias. First 
the Hubbard interaction is switched off for the test of 
the calculation precision for a chain of L = 40 without 
the special bond. The lengths of ideal leads attached the 
chain are taken to be long enough so that the results are 
sufficiently insensitive to it. The incoming and outgoing 
currents are defined as the corresponding currents at the 
left and right interfaces according to Eq. (3), respectively. 
At U = 0, the spin and charge currents are the same since 
no interaction between electrons of different spins exists. 
Furthermore, the currents can be calculated exactly at 
the case, which is shown as solid lines in (a). A time lag 
between the incoming and outgoing currents is observed 
due to the spin/charge transportation from the left chain 
end to the right one. Finally, steady transport is reached 
at t = i0(h/to). In Fig. 4(a), we also show the results 
obtained from the t-DMRG method by keeping various 
numbers of DMRG states. It can be seen that the curves 
obtained when the number of kept DMRG states (M) is 
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FIG. 5: The dependence of transmission rates (see text for 
definition) on (a) tj, (U = St ) and (b) U {t\ = 0). 



128 or more are very close to the exact one (solid lines), 
while the convergence for the outgoing current is faster 
than that of the incoming one, the reason being that the 
voltage bias is applied only at the left interface for con- 
venience in this calculation. 

Next, in Fig. 4(b) we show both spin and charge outgo- 
ing currents for a Hubbard chain of L = 40 and U = 8io 
with the special bond for a number of different kept 
DMRG states (M). It is clearly seen that the conver- 
gence for a Hubbard chain is much better than that of the 
U = chain with increasing the number of kept DMRG 
states (M). In a sufficiently long time, it can be seen 
that while the charge current passes through the chain 
the spin current is blocked and only a little of it passes. 
This result is consistent with the spread of a wave packet 
studied above. In Fig. 4(c) we give the currents for dif- 
ferent lengths (L). We observe that a Hubbard chain of 
moderate length L (= 30 — 40) is optimal, for a shorter 
chain is accompanied by a serious finite-size effect and a 
longer chain by a large accumulated error since it takes 
more time for spin/charge to travel from one end to the 
other. 

Finally, we give the dependence of spin/charge trans- 
port through the Hubbard chain on the parameters t± 
and U in Fig. 5. The transmission rates are defined 
as the ratio between transmitted and incident currents, 
where the former is taken from the average value of the 
outgoing current at a period when the current is almost 
steady and the latter is the same of the incoming current 
at the period before the current reaches the right end. 
From Fig. EJa), it can be seen that the transmission rate 
changes little for t± < 0.5. On the other hand, Fig. [5Jb) 
shows that the on-site Hubbard repulsion U plays a cru- 
cial role in the observed spin-current blockade. When U 
increases, the charge transmission rate changes little in 
the strongly-correlated regime, namely, U > 4irj, and the 
spin rate decreases very quickly and reaches about 0.1 for 
U > IUq. It implies the spin current is nearly completely 
blocked by the special bond. 

An experimental observation of spin-charge separation 
in cold Fermi gases has been proposed by Kollath, ct 
al.,[23[ in which the energy scale fcsT should be much 
smaller than the Mott energy gap to ensure that thermal 



activation does not destroy the Mott-insulating behavior. 
But this is not required in our case since the Mott gap 
is zero. Furthermore, an important advantage is that 
the spin-current blockade could be realized by adjusting 
only one parameter (fj,), which simplifies experimental 
observation. 

In summary, we have investigated the spin and charge 
transport in a Hubbard chain with a bond of control- 
lable electronic hopping. We find the spin current can be 
blocked by this special bond while charge current passes 
through the bond freely. It is found that a large Hub- 
bard U is required for the observed blockade since the 
spin-charge separation plays a crucial role in it. 
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